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Abstract 

It is shown that the fusion of strings is a source of particle production in 
nucleus-nucleus collisions outside the kinematical limits of nucleon-nucleon 
collisions. The spectrum of different particles is compared with the high 
energy data on p-A collisions obtaining a reasonable agreement. Results 
for A-B collisions at yfs = 19.4 AGeV and ^fs = 200 AGeV are given. It is 
shown that the fusion of strings can accelerate particles up to the highest 
energy detected in cosmic rays without help of any additional cosmic ac- 
celerator. Also the rise of the average shower depth of maximun for cosmic 
rays in the energy range between 10 16 eV and 10 19 eV can be explained by 
the same mechanism without requiring any change in the chemical compo- 
sition. 
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1 Introduction 



Several models of hadronic interactions have been very successful in describing particle 
production in hadron-hadron, hadron-nucleus and nucleus-nucleus collisions. Monte 
Carlo versions of these models are in reasonable agreement with most of the prop- 
erties of soft multiparticle production. In these models (|], |2], |3|, [|]) strings, chains 
or pomerons are exchanged between the projectile and target. The number of strings 
grows with the energy and with the number of nucleons of the participant nuclei. In the 
first approximation strings fragment into particles and resonances in an independent 
way. The only correlation among strings is due to energy-momentum conservation. 
However, the interaction between strings becomes important with their number grow- 
ing. Even at SPS energies both a large number and a large density of strings are 
expected. For instance, for S-S and Pb-Pb central collisions the estimated number 
of strings are 120 and 1300, respectively, and their densities 3.3 fm~ 2 and 9.8 fm~ 2 , 
respectively. At the relativistic heavy ion collider (RHIC) and large hadron collider 
(LHC) energies the number and density of strings are not negligibly small even for 
hadron-hadron collisions. 

The interaction between strings or the interaction of resonances produced in the 
fragmentation of strings have been introduced in some of the models ( H |6[, |7| |8|, |9j, pT0f| ) . 
In particular, fusion of strings has been incorporated into the Dual Parton Model 
(DPM) (U) and the Quark Gluon String Model (QGSM) ((JTTJ). Some of the effects 
of string fusion, like strangeness and antibaryon enhancement ([12|]), reduction of long 
range correlations (|13[) an d multiplicity suppression, are in reasonable agreement with 
the existing experimental data. Also predictions for the Relativistic Heavy Ion Collider 
(RHIC) and Large Hadron Collider (LHC) are avalaible. In this paper we explore 
another effect of string fusion, namely particle production in nucleus-nucleus collisions 
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outside the kinematical limits of nucleon-nucleon collisions (the so-called cumulative 
effect). 

It is shown that at present avalaible energies a non negligible number of baryons 
and mesons are produced with momenta greater than the ones of the colliding nucle- 
ons. Essentially, the particles produced outside the kinematical limit are protons and 
neutrons but an appreciable number of 7r's and A's are also predicted. This effect, 
together with the reduction of multiplicities, provides a natural explanation of some 
features of cosmic ray data, like the rise of the average shower depth of maximum 
X max (the amount of air penetrated by the cascade when it reaches maximum size) 
(EH 15 1 ) with increasing energy from 10 17 eV to 10 19 eV, and the existence of events 



with energy above 10 20 eV (@, [R| jT^j), higher than the expected cut-off (JT|, |9|, |0 



due to the scattering of cosmic rays with the microwave radiation background. Usu- 
ally the first feature is explained by an enrichment of protons in the composition of 
primary cosmic rays ([|15|]) as energy increases. However, as we shall show, if the com- 
position of the primary cosmic rays is kept fixed in the energy range between 10 17 eV 
and 10 19 eV, string fusion leads to a suppression of the multiplicity similar to the one 
produced by changing heavy nuclei (Fe) by protons in the composition of the primary. 
On the other hand, since the momentum of a fused string is a sum of the momenta of 
its ancestor strings, it is possible to obtain particles with more energy than the initial 
nucleon-nucleon energy. As in a considerable part of events the primary would be iron, 
several string could fuse in central Fe-Air collisions. Therefore, the observed cosmic 
ray events with energy above 10 20 eV may actually correspond to three or four times 
less initial energy per nucleon than the one apparently measured. Reducing the energy 
by a factor three, the attenuation mean free path can grow more than an order of 
magnitude. String fusion could make these events compatible with the existence of the 
above mentioned cut-off. 



2 String fusion in the Monte Carlo code 

To study the particle production outside the kinematical limits of nucleon-nucleon 
collisions we use a Monte Carlo code based on the QGSM, in which the fusion of 
strings has been incorporated (fll|). A detailed description of the Monte Carlo String 
Fusion Model (SFMC) and comparison with experimental data can be found in Refs. 



TT| , 12, 13 1 . A hadron or nucleus collision is assumed to be an interaction between 
clouds of partons formed long before the collision. Without string fusion partons are 
assumed to interact only once. Each parton-parton interaction leads to the creation of 
colour strings. Since both the projectile and the target must remain colourless, strings 
have to be formed in pairs. For instance, in nucleon-nucleon collisions at moderate 
energies a pair of strings is formed between a valence diquark of the projectile and a 
valence quark of the target and between a valence quark of the projectile and a valence 
diquark of the target. As the energy increases, pairs of strings are also formed between 
sea quarks (antiquarks) of the projectile and sea antiquarks (quarks) of the target. 
Hadrons and nuclei are considered on the same footing. The nuclear wave function is 
taken as a convolution of the parton distribution in a nucleon with the distribution of 
nucleons in the nucleus. For the distributions of partons and nucleons we take gaussian 
(centered in each nucleon) and Wood-Saxon shapes respectively. 

Strings fuse when their transverse positions come within a certain interaction area, 
which is fixed previously to describe correctly the strangeness enhancement ([0) shown 
by the data on nucleus-nucleus collisions. Fusion can take place only when the rapidity 
intervals of the strings overlap. It is formally described by allowing partons to interact 
several times, the number of interactions being the same for projectile and target. The 
quantum numbers of the fused string are determined by those of the interacting partons 
and its energy-momentum is the sum of the energy-momentum of the ancestor strings. 
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The colour charges of the fusing string ends sum into the colour charge of the resulting 
string ends according to the SU (3) composition laws. In particular, two triplet strings 
fuse into an antitriplet and a sextet string, with probabilities 1/3 and 2/3 respectively. 
A triplet and an antitriplet string give rise to a singlet state and an octet string with 
probabilities 1/9 and 8/9 respectively. In present calculations only fusion of two strings 
is taken into account. 

A quasi-classical picture of the decay of colour strings is assumed in which pairs of 
oppositely colour partons are produced in the string colour field, which neutralizes this 
field and leads to string breaking. The new sextet and octet strings are supposed to 
break with the production of two (anti) quark complexes with the same colour charges 
Q and — Q as those of the ends of the string. The created (anti) quarks have arbitrary 
flavours and masses chosen as the corresponding constituent masses. The probability 
rate for the constant colour field of two opposite charges Q-Q to create a parton pair 
with the same colour charges Q-Q and transverse mass M t for unit string length and 
time is taken by the Schwinger expression 

W ~ K* N] exp(-nM?/K [N] ) , (1) 

where K[ N ] is the string tension for the [N] SU (3) representation proportional to the 
corresponding quadratic Casimir operator C^ N y In our case 

Cf 3] = 4/3, Cf 6] = 10/3, Cf 8] = 3 . (2) 

The fragmentation of the fused strings produces more baryons and antibaryons, es- 
pecially strange ones, than in the case of the fragmentation of the original strings. 
These enhancements are accompanied by a strong reduction of particle production in 
the central rapidity region. 
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3 Cumulative particle production 

Particle production outside the nucleon-nucleon kinematical limits is a well known 
effect, called cumulative effect, studied both theoretically and experimentally ([^, [23|, 



2% |25| , p6|). However at high enough energies, where the string picture can be applied, 
there are only data from one collaboration at 400 GeV/c (|[26||), with incoming protons 
against nuclei: Li, Be, C, Al, Cu and Ta. We have generated 10000 events for each p-A 
collision in our Monte Carlo code. The results together with the experimental data for 
the invariant differential cross sections for production of protons, positive pions and 
positive kaons are shown in Tables 1 and 2. 

If the invariant differential cross section is parametrized as a function of xp in the 
exponential form: 

a ~ exp(-bx) (3) 

the slope b is in the range 5-6 for all cases in agreement with the experimental data. 
If the fusion of strings is not included the obtained b value is in the range 10-11 quite 
far from the data. 

The fusion of strings describes also rightly the dependence of the inclusive cross 
section with the kinetic energy. As an example, in Fig. 1, it is shown the experimental 
data for the positive pion inclusive cross section for p-Ta collisions together the results 
with and without fusion of strings. 

The A-dependence of the cumulative effect is usually parametrized as A a ^ x \ The 



value of a(x) has been proposed (p4j) to be the function 

(x~A Q(x) , a{x) = 1 + 1/3 (x-1). (4) 

which is plotted in Fig. 2 together the a(x) obtained from the experimental data which 
is very close to the above straigth line. Also the a(x) obtained from the fusion of strings 
is close to the experimental data, indeed it fits into the line a(x) = 1 + 0.4(x — 1). 



From all these comparisons we can conclude that a reasonable agreement for p-A 
cumulative particle spectrum is obtained. Notice that we do not have any free param- 
eter. We could improve the comparison by introducing minor effects like rescattering 
or a more detailed nucleon and quark correlations (other than Fermi motion, which 
has already been included in our code). However our goal is not to obtain a perfect fit 
but just to check if the string fusion works reasonable well. 

Once it has been checked that our model is consistent with the hadron-nucleus 
experimental data, we turn to nucleus-nucleus collisions , simulating 10000 central S- 
S collisions and 1000 central Pb-Pb collisions at yfs — 19.4 AGeV. Also central Pb-Pb 
collisions at RHIC energies (yfs = 200 AGeV) have been simulated. Distributions of 
baryons and mesons in central S-S collisions at SPS energies with xf larger than 1 are 
shown in Fig. 3. In Fig. 4 the separated spectra of protons, neutrons, lambdas, K°, 
K + , K~ , 7T°, 7r + and n~ are presented. Fig. 5 and Fig. 6 show the same distributions 
for the case of central Pb-Pb collisions at yfs = 19.4 AGeV. The results for Pb-Pb 
collisions at = 200 AGeV are very similar to the ones at yfs = 19.4 AGeV. In 
1000 events 2015 particles are found with \xf\ larger than 1 to compare with 1783 at 
yfs = 19.4 AGeV. This small change is due to a moderate increase of the number of 
strings with energy. In general, it is seen that many baryons and mesons are produced 
with xf larger than 1. 



4 String fusion at ultrahigh cosmic ray energies 

String fusion produces a strong supression of multiplicities. In the limit of a very strong 
fusion, the multiplicity in hadron-nucleus collisions turns out independent of A instead 
of ~ A 1 / 3 . Also the squared dispersion D 2 behaves like A' 1 / 3 instead of A 1 / 3 . At finite 
energies the fusion of strings is not strong and it is not expected such a behaviour, 



however the reduction of multiplicities is quite sizable for heavy ion collisions at RHIC 
and LHC energies. For instance, simulations done for ALICE detector for LHC indicate 
that the multiplicity is a factor between 2.5 and 4 less than the one obtained in the 
models without fusion. 

This suppression of multiplicities can explain the rise of the average shower depth 
of maximum in cosmic rays as the energy increases, without requiring any change in 
the chemical composition. It is usually accepted that there is a change in the cosmic 
ray chemical composition between 10 16 eV and 10 19 eV. It seems that the composition 
becomes significantly lighter with increasing energy, going from a heavy composition at 
10 16 eV to a light one at energies higher than 10 19 eV. The distribution of the shower 
depth of maximum as a function of energy has been studied using a simple model of 
two components ( ||15|| ), observing that the composition of the primary changes from 
approximatly 75 % of iron component and a 25 % of proton component at 10 16 eV to 
50 % of iron and a 50 % of proton at 10 19 eV. To study this point, we have computed 
the multiplicities of p-Air and Fe-Air interactions with and without string fusion in 
the whole range of energies studied (from 10 16 to 10 19 eV). As it can be seen in Fig. 7, 
with string fusion the multiplicity for a constant composition of 10 % of proton and 90 
% of iron in the whole range of energy, essentially reproduces the multiplicity obtained 
without string fusion for a uniform change in the composition from 75 % Fe and 25 % 
proton at 10 16 eV to 50 % Fe and 50 % proton at 10 19 eV. Thus the string fusion does 
the same job as the composition change. 

Therefore, the change in the energy behaviour of the average shower depth of maxi- 
mum X max can be ascribed to a change in the interaction mechanism with the increasing 
role of collective effects like string fusion, and not to a change in the chemical composi- 
tion of the primary cosmic rays. Further studies of this point would require combining 
the code used in this paper with the standard codes which describe the full cascade. 
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Work in this direction is in progress. 

Once we know that it is reasonable to assume that most of the primary cosmic 
rays are iron nuclei, we would like to know how many particles are going to be accel- 
erated due to the string fusion mechanism and get energy-momentum larger than the 
permitted limits in nucleon-nucleon collisions. 

To study the case relevant for cosmic rays we simulated 1000 Fe-Air collisions at 
10 17 eV (we used this energy and not 10 20 eV to save computing time, rendering the 
simulation reliable). In this sample 198 particles with \xp\ > 1 were found. The average 
number of strings was found to be 225, from which 62 joined to form double strings. 
As mentioned, our code only includes fusion of two strings. However we can estimate 
the number of strings participating in a triple fusion assuming that the probability for 
triple fusion is roughly the square of that for double fusion. Then one would expect 
that 18 strings join to form triple strings and 4 strings join to form a quadruple string. 
Other reasonable assumptions about the probability of triple fusion give similar results. 
Therefore the probability of obtaining particles with \xp\ > 2 or even \xf\ > 3 does not 
seem to be negligible. Triggering central collisions, particles with \xp\ > 4 could even 
be detected. The energy around 3 • 10 20 eV measured in several cosmic ray experiments 
could then be lowered by a factor 2 to 4 if the described effect is present and there 
are particles in the shower with \xp\ > 2 or \xp\ > 3. This lower energy for the 
primary may lie below the energy cut-off due to the scattering of cosmic rays on the 
microwave background. Indeed, according to the computation of F. W. Stecker d|T9|) 
the attenuation mean free path increases more than one order of magnitude. This fact 
means that the measured energies are compatible with the effective cut-off. 
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5 Conclusions 



The string fusion mechanism produces particles with more energy-momentum than 
the original nucleon-nucleon collisions. In this way it can be considered as a hadronic 
accelerator. We have shown that string fusion reproduces reasonably well the experi- 
mental data on cumulative effect in proton-nucleus collisions giving a sizable number 
of protons, neutrons, lambdas and pions in central nucleus-nucleus collisions already 
at t/s = 19.4 AGeV. The experimental confirmation of this at SPS will be welcome. 

Concerning cosmic rays, string fusion can explain the reduction of multiplicities in 
the energy range between 10 16 eV and 10 19 eV without requiring a change in the chem- 
ical composition. Also, string fusion provides a natural hadronic accelerator to reach 
energies above 10 20 eV without requiring any additional unusual cosmic accelerator. 
Our predictions can be checked in future heavy ion experiments at the accelerators 
RHIC, LHC and also in cosmic ray experiments (concretely the Auger proyect (P7[)). 

In conclusion we would like to thank N. S. Amelin, A. Capella, J. W. Cronin, G. 
Parente, J. Ranft and E. Zas for useful comments and discussions and the Comision 
Interministerial de Cienci'a y Tecnologia (CICYT) of Spain for financial support. M. 
A. Braun thanks IBERDROLA, E. G. Ferreiro thanks Xunta de Galicia and Yu. M. 
Shabelski the Direction General de Pohtica Cientifica of Spain for finantial support. 
This work was partially supported by the INTAS grant N 93-0079. 
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Table captions 



Table 1. Comparison of experimental data ( ||26|| ) on the invariant differential cross 
section a = E^ (GeV mb/(GeV/c) 3 sr nucleon) for n + and K + vs p (GeV/c), labo- 
ratory angle 118°, P[ a b = 400 GeV, for p-Li and p-Ta collisions with the String Fusion 
Model code results, with and without string fusion. 

Table 2. Comparison of experimental data ( ||26|| ) on the invariant differential cross 
section a = E-^ (GeV mb/(GeV/c) 3 sr nucleon) for n + vs p (GeV/c), laboratory 
angle 118°, P[ a b = 400 GeV, for p-Li and p-Ta collisions with the String Fusion Model 
code results, with and without string fusion. 
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Figure captions 



Fig. 1. Comparison of experimental data ([26]) (continuous line) on the invariant 
differential cross section a = for n + vs kinetic energy (GeV), laboratory angle 

118°, Pi a b = 400 GeV, for p-Ta collisions with the String Fusion Model code results, 
with (upper dashed line) and without string fusion. 

Fig. 2. a vs x. a = 1 + 1/3 (x — 1) (continuos line), experimental result (nearer 
dashed line) and our results with fusion. 

Fig. 3. xp distributions for xf > 1 in S-S collisions (10000 events) at y/s = 19.4 
AGeV of mesons (a) and baryons (b) with (continuous line) and without (dashed line) 
string fusion. No mesons are found in the no fusion case. 

Fig. 4. xf distributions for xp > 1 in S-S collisions (10000 events) at y/s = 19.4 
AGeV of protons (a), neutrons (b), lambdas (c), K° (d), K+ (e), K~ (f), tt (g), tt+ (h) 
and ix~ (i) with (continuous line) and without (dashed line) string fusion. No mesons 
are found in the no fusion case. 

Fig. 5. xp distributions for xp > 1 in Pb-Pb collisions (1000 events) at y/s = 19.4 
AGeV of mesons (a) and baryons (b) with (continuous line) and without (dashed line) 
string fusion. No mesons are found in the no fusion case. 

Fig. 6. xp distributions for xf > 1 in Pb-Pb collisions (1000 events) at y/s = 19.4 
AGeV of protons (a), neutrons (b), lambdas (c), K° (d), K + (e), K~ (f), 7r° (g), rt + (h) 
and 7T~ (i) with (continuous line) and without (dashed line) string fusion. No mesons 
are found in the no fusion case. 

Fig. 7. Total multiplicity dependence on the primary energy for a fixed composition 

< n t > = 0.1 < n p _Air > + 0.9 < n Fe -Air > in the fusion case and a uniform change 
in the composition from < n t > = 0.25 < n p _Air > + 0.75 < np e -Air > at 10 16 eV to 

< n t > = 0.5 < n p _Air > + 0.5 < n Fe _Air > at 10 19 eV in the no fusion case (dashed 
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line). 
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Table 2 
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